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Preliminary Sunyaev Zel'dovich Observations of Galaxy 
Clusters with OCRA-p 
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ABSTRACT 

We present 30 GHz Sunyaev Zel'dovich (SZ) observations of a sample of four galaxy 
clusters with a prototype of the One Centimetre Receiver Array (OCRA-p) which 
is mounted on the Torun 32-m telescope. The clusters (CI 0016+16, MS 0451.6-0305, 
MS 1054.4-0321 and Abell 2218) are popular SZ targets and serve as commissioning 
observations. All four are detected with clear significance (4 — 6a) and values for the 
central temperature decrement are in good agreement with measurements reported 
in the literature. We believe that systematic effects are successfully suppressed by 
our observing strategy. The relatively short integration times required to obtain these 
results demonstrate the power of OCRA-p and its successors for future SZ studies. 

Key words: cosmology: observations - cosmic microwave background - galax- 
ies:clusters:individual (CI 0016+16, MS 0451.6-0305, MS 1054.4-0321, A2218) - meth- 
ods:observational 



1 INTRODUCTION 

The SZ effect is a spectral distortion of the Cosmic Mi- 
crowave Background (CMB) caused by the inverse Comp- 
ton scattering of CMB photons off the hot plasma found in 
clusters of galaxies. At low radio frequencies the SZ effect 
manifests itself as a fractional decrement in the CMB of or- 
der 10 -4 , whereas at frequencies greater than ~ 220 GHz an 
increment is observed. 

Since the f irst SZ detections in 
l|Birkinshaw et all Il978l ) much progress has 
Extensive interferometric s tudies have been 
from the first SZ images l|jones et alj 1 19931 ). to the use 
of cluster sample s to perform cosmological studies (e.g 



the 1970s 
been made, 
undertaken, 



20011 . iReese et all |2002|. ISaunders et al l 12001 
20051 and | Bonamente et al.ll2006l ). Radiometers 
(e.g. Myers et al.lll997l ) and bolometers fe.g. iBenson et all 



Greg o et al 
Jones et alj 



reviews of lBirkinshav^lll999l and ICarlstrom et al.ll2002T ). Al- 
though most work in the field of SZ studies to date has been 
performed using non-ideal instruments, numerous purpose- 
built SZ observatories ar e current ly under constructio n (e.g. 
AMiBA JLo et all [200ch AMI feneissl et all l200ll ) SZA 
llLoh et all 120051 ). ACT (|Kosowskvl 120031 ) and SPT (|Ruhll 
12004 ) ) all with the common aim of performing blind surveys. 
Such work will exploit the unique redshift-independence of 
the SZ effect to produce catalogues which are essentially 
mass-limited, and thus less affected by the biases which 
plague selection via optical or X-ray methods. Studies of 
number counts will be used to further constrain the cosmo- 
logical model and study cluster evolution from formation 
until the present. Those telescopes which have the capabil- 
ity to produce high-quality imaging will additionally provide 
the opportunity to better understand cluster physics. 



|2004| ) have also been used to good effect. More recently, 
purpose-built C MB instruments have proved their worth 
in S Z studies (jLancaster et alj 120051 . lUdomprasert et all 
|2004| ). although contamination from primordial anisotropies 
can severely limit the achievable signal-to-noise ratio. 

Measurements of the thermal and kinematic SZ effects 
are studied for the unique information that they can provide 
on cosmology and the structures of cluster atmospheres (see 



OCRA (|Browne et all |2000| ) is a planned 100-element 
continuum receiver which will have excellent surveying and 
imaging capabilities, and will thus be ideal for performing 
blind surveys in order to study radio source populations 
and SZ clusters. The prototype for OCRA, OCRA-p, is a 
two-element receiver mounted on the 32-m telescope at the 
Torun Centre for Astrophysics of the Nicolas Copernicus 
University, Poland. Even at this preliminary stage, while the 
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instrument concept is being tested, high-sensitivity SZ mea- 
surements are possible, with imaging and surveying planned 
for the future. We here present our first observations of 4 
well-studied massive SZ clusters (CI 0016+16, MS 0451.6- 
0305, MS 1054.4-0321 and A2218). 

The structure of this paper is as follows. In Section 2 we 
give a brief description of the Torun 32-m telescope and the 
OCRA-p receiver. Section 3 contains details of the cluster 
sample. The observations and data reduction are described 
in Section 4, and the problem of radio source contamina- 
tion is detailed in Section 5. Section 6 presents our results 
and discussion. Throughout the paper we adopt the follow- 
ing cosmological parameter values: Ho = 70kms~ 1 Mpc" 1 , 
fi m0 = 0.3, fi A o = 0.7 



2 THE TORUN TELESCOPE AND OCRA-P 

The Torun observatory is located in Piwnice, 15 km out- 
side Torun in northern Poland. The telescope consists of a 
32-m parabolic dish and 3.2-m sub-reflector, with a fully 
steerable classical Alt-Az mount. It has receivers operat- 
ing at 1.4-1.7, 5, 6.8, 12 and now, with OCRA-p, 30 GHz. 
The telescope is used for a variety of studies incl uding inter- 



stell ar molecules (e.g. Blaszkiewicz fc Kusll2004l ). and VLBI 



(e.g. iBartkiewicz et alj|2005l ) 

The OCRA prototype, OCRA-p, is funded by a grant 
from the Royal Society Paul Instrument f und. The instru- 
ment is described in detail in lLowei l|2005h and lLowe et al.l 
(|2007h ; here we present a short summary. The basic radiome- 
ter design is based on the prototype for the Planc k Low 
Frequency Instrument (LFI, Mandolcs i et alj (120001)), and 
is sim ilar to the WMAP K-band receivers l|Jarosik et al.l 
2003). OCRA-p consists of two horn-feeds, the beams of 
which are separated by 3'.1 and have FWHM 1'.2. As the 
beam separation is small, it is possible to reduce the effects 
of atmospheric and gain fluctuations by switching between 
the input channels and taking the difference. This can be im- 
proved upon by further levels of switching. The full switching 
strategy for SZ observations is described in Section \4. II 



3 THE CLUSTER SAMPLE 

For these preliminary observations, we have chosen four clus- 
ters which have well-studied SZ effects, and which cover 
a wide range of declinations, in order to assess the per- 
formance of the telescope/receiver combination. In a fu- 
ture paper we will present results from observations of a 
large, statistically-complete sample of clusters, observations 
of which are underway. The cluster sample, and its X-ray 
properties, is summarised in Table [T] 

Our four chosen clusters have been studied at radio 
wavelengths by a variety of groups employing wide-ranging 
techniques. In the subsequent sections we summarise pre- 
vious work in terms of the most recent, any that is close 
to our observing frequency, and any that employs a similar 
observing strategy. These clusters have also been studied in 
detail in X-rays. We note briefly the interesting character- 
istics reported in the literature, paying particular attention 
to those which may affect the SZ decrement. We also report 



detections of lensing by our target clusters where appropri- 
ate. 



3.1 CI 0016+16 

CI 0016+16 was observed by the pioneers of SZ astron- 
omy, and continues to be a popular target. The SuZIE 
team performed simultaneous bolometric observations in 
three frequency bands (145, 221 and 355 GHz) and find a 
central Comptonisa tion parameter j/o of (3.2712'gei) x 4 



dBenson et aUl2004h. t hus AT = -1785t™i 6 MK- The work 
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Tsuboi et al.1 (|2004h is interesting, as they employ sim- 
ilar observing techniques to those described in this work. 
They find a central decrement of —1431 ± 133 /iK us i ng the 
Nobeyama 45-m telescope, at 43 GHz. iReese et all (|2002l ) 
find a similar result of —1242 ± 105 /iK, observing with the 
OV RO and BIMA interferometers at 30 GHz. 

IWorrall"fc~ Birkinshaw (2003) present analysis of an 
XMM-Newton observation of CI 0016+16. They find a best 
fitting temperature of 9.131q '22 keV within 1.5arcmin of 
the cluster centre (see Table [TJ, an d note the existence 
of a central structure which may have a harder spec- 
trum. This, along with the presence of an emission re- 
gion to the west, may be evidence of merger activity. They 
comment that their determination of the total gravitating 
mass is in g ood agreement wi th the value derived from 
lensing data l|Smail et al.l I1997T ). thus supporting assump- 
tions of hydrostatic equilibrium. Combination of their X- 
ray results with a reanalysi s of the SZ data presented in 
iHughes fc Birkinshawl (|l998l ) yields a value for the Hubble 
constant of 68 ± 8kms _1 Mpc _1 (random error only). 



3.2 MS 0451.6-0305 



iBenson et all j|2004l ) also observed MS 0451.6-0305, find- 
ing a central Comptonisation parameter from their multi- 
frequency data of (2.84 ± 0.52) x 10 -4 , corresponding to 
ATq = —1550 ± 284 /jK. This c luster was also included in 
the recent work o f Tsub oi et al. I l|2004l) . who record a central 
decrement of -1201 ± 184 ^K. MS 0451.6-0305 also features 
in the work of the OVRO /BIMA collaboration. IReese et al.l 
(2002) quote a central decrement from their 30 GHz mea- 
surements of — 1431!tg3 ^K. 

The gas temper ature and luminosity o f MS 0451.6-0305 
are given in Table [T] iDonahue et al.l (|2003l ) analyse a Chan- 
dra observation in order to obtain a single temperature. 
They extract spectra from annular regions and do not find 
a significant gradient between the inner and outer regions 
of the cluster using these data. However, they also find an 
acceptable fit to a two-temperature plasma model, although 
the associated cool component would not c ontrib ute signif- 
icantly to the X-ray emission. ISand et al.l (|2005l ) identify 6 
candidate lensing arcs in MS 0451.6-0305. 



3.3 MS 1054.4-0321 

MS 1054.4-0321 is less w ell-studied at radio w avelengths. It 
was observed recently by [Benson et al.l (|2004l ). who find the 
central Comptonisation parameter to be yo — i^-^-i'it) x 
10~ 4 from the measurements at 145 GHz. This corresponds 
to a temperature decrement of 2113lgi2 mK- This cluster 
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Cluster 


RA 


DEC 


z 


Lx 


T x 


References 




(J2000) 


(J2000) 




(lO^ergs" 1 ) 


(keV) 




CI 0016+16 


00 18 33.30 


+16 26 36.00 


0.546 


22.7 


9 13 +UM 


j, h, i 


MS 0451.6-0305 


04 54 10.90 


-03 01 07.00 


0.550 


17.5 


8.05l°i 


b, f, g 


MS 1054.4-0321 


10 57 00.20 


-03 37 27.00 


0.823 


12.3 


7 O+0.7 
'•^-0.6 


b, e, e 


A2218 


16 35 49.10 


+66 12 30.00 


0.176 


9.6 


6 63+ ' 27 


a, c, d 



Table 1. Basic cluster data. Coordinates correspond to pointing centres, para meters are as listed in the BAX X-ray cluster data base. 
Luminosities refer to the 0.1-2.4 keV R OSA T energy band. Refe rences: (a) | Struble fc Rood! Il999j. (b) iGioia fc Luppind Il994l, (c) 
iBalogh et alll2002i. (d) iPratt et al.ll2005l. (e) iGioia et al.||2004 (f) iDonahue et al.ll2003l. (el iTozzi et al.ll2003l (hi lEllis fc Jonesl 12003. 
(il lWorrall fe Birkinshawll2003l , (i) lDressler fc Gunnlll992l 



has als o been obs e rved at 30 GHz using the OVRO/BIMA 
arrays. IJov et all (|200lf ) use the SZ d ata to infer a clus - 
ter temperature of 10.4l^o keV, whereas [Grego et al.l (120011 ) 
calculate the gas mass fraction within rsoo (where rsoo is the 
radius inside which the cluster density is greater than 500 
times the critical density) to be / g = 0.053 ± 0.028. 

iJee et all (|2005T ) present a joint analysis of Chandra and 
HST weak lensing data for this cluster. They find three dom- 
inant clumps in the dark matter distribution, and observe 
that one substructure is not present in the X-ray data. The 
remaining two features are displaced in the X-rays relative 
to the HST shear map, probably as a result of a current 
merger. They find consistent mass estimates (r < IMpc), 
using the datasets independently, of ~ 1 x 10 Mq . They 
also report a significant temperature gradient, but comment 
that no shock-heated region is observed between the X-ray 
peaks. 



3.4 Abell 2218 

A2218 was one of the first clusters to be detected via the SZ 
effect in the 19 70s, and has been stu died in much detail since. 
Most recently, IJones et all |2005) observed the cluster at 
15 GHz with the Ryle telescope, and found a maxim um cen- 
tral decrement of -760±150 £iK.|R eese et all (|2002 J ) reco rd a 
central decrement of -73lt {H uK. iTsuboi et all (Il998l ) ob- 
served this cluster, this time at 36 GHz, and measured the 
SZ effect to be -680 ± 190 /iK. 

A2218 has been studied extensively at X-ray wave- 
lengths, and has an X-ray luminosity of 9.55 x 10 44 ergs -2 
and an average gas temperature of 6.63 ± 0.27 keV. A Chan- 
dr a study of BCS c lusters including A2218 is presented 
in iBauer et al.l l|2005h . They find that the cluster exhibits 
a cool core. This cluste r is also part of a HST lensing 
study l|Smith et al.ll2005[ ). The authors combine strong and 
weak lensing information in order to deduce a total cluster 
mass of (5.6 ± 0.1) x 10 14 Mq. Using additional information 
from archival Chandra observations, they compare the mass 
and X-ray morphology of the cluster to establish quantita- 
tively that it appears dynamically immature. This is sup- 
ported by other Ch andra studies, such as that reported in 
iGovoni et al.l f2004). This work reinforces previous claims 
that A2218 has an irregular temperature structure, despite 
its relatively symmetric X-ray morphology. 



4 OBSERVATIONS AND DATA REDUCTION 
4.1 Observing Strategy 

The observations presented in this paper were performed at 
intervals between June 2005 and April 2006. We aimed to 
achieve a similar noise level for all clusters, thus observing 
times were roughly 11 hours for each. The time spent on the 
cluster itself was less than half of this value due to telescope 
control issues which have since been resolved. Changes in 
the atmosphere are problematic for single-dish experiments, 
and we employ a differencing technique in order to remove 
the resulting effects. For each cluster we make two observa- 
tions, firstly the cluster itself (hereafter referred to as the 
TARGET field) and secondly a trailing field (hereafter re- 
ferred to as the TRAIL field) at the same declination but 
offset in right ascension. To this TARGET-TRAIL strategy 
we apply both beam switching and position switching. 

Firstly for the TARGET field, the horn-feeds are po- 
sitioned such that one beam, beam A, is coincident with 
the cluster centre and the other, beam B provides a mea- 
sure of the blank sky signal. (For extended sources, beam 
B may measure a small signal itself. This must be properly 
accounted for - see Section [4. 41 ) We switch between the A 
and B beams at a rate of 280 kHz, recording the integrated 
A — B difference every second. The beams are repositioned 
every 53 seconds such that beam B measures the cluster 
and beam A measures the sky background, and the differ- 
encing is repeated. We then take the double-difference, i.e. 
(A\ —Bi) — (£?2 —A2), such that we recover twice the cluster 
signal relative to the two background regions. This process 
is then repeated for a TRAIL field, which is observed over 
the same range of hour angle as the TARGET field. We 
are satisfied that our TRAIL fields contain at worst negligi- 
ble amounts of contaminant signal and thus retain the data 
purely as a check. 



4.2 Calibration 

The data are calibrated against an internal noise source, 
which is itself calibrated via observations of the well-known 
bright radio sources NGC7027 and 3C286 with assumed flux 
densities 5.64 Jy El and 2.51 Jy. Significant changes in the 
level of the internal noise source occurred over the dura- 
tion of these observations. Some were due to changes in the 



1 We note tha t this value, extrapolated to 30 GHz via the scale 
of iBaar s et alj lll977T ). is 3.5% lower than the measurements of 
iMason et al. I 1l999)) . This is well below our systematic error level 
of 5%. 
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way that the noise source was used, but others were due 
to a temperature sensitivity of the noise source, which was 
controlled by improving the temperature stability in the re- 
ceiver cabin. Residual uncertainties in the calibration of the 
system at the level of 5% may be expected. 

4.3 Statistical Data Analysis 

After combination of the second-by-second average data into 
double-differenced measurements of the brightness of the 
sky, and calibration, the data are examined for periods of 
increased noise (which might arise from receiver instabilities 
or bad weather conditions) or individual anomalous points. 
The quality of the data throughout the observing period was 
good, with a tendency to show a slight increase in noise in 
Spring 2006 as the weather became warm and wet: weather 
effects were minimized by the bulk of the observations being 
taken in visually clear conditions. 

The combination of the data into final averages was 
performed including statistical tests for outlier data. The 
fractions of the data points rejected by 3a or 5a cuts were 
small in all cases, and no cut-dependent changes in the aver- 
age results were seen. The distributions of data values in the 
double-difference data are close to Gaussian, with a slight 
tendency to show elevated wings in the distributions: the 
estimates of the error on the mean (Table [4) take account 
of these deviations from Gaussian distributions. 



4.4 Removal of Instrumental Signature 



As mentioned in Section 14.11 the small separation of the 
OCRA-p beams relative to the extent of a cluster's SZ decre- 
ment leads to a significant reduction in the measured signal 
compared to the intrinsic central signal. In order to remove 
this effect and recover the 'true' SZ signal, we adopt the 
standard /3-model approach. We take values for the core 
radius r c and f3 from the literature (see Table [2} in order 
to estimate the fraction, F of the SZ signal measured by 
beams A and B respectively, given where the beams fall 
on the sky, and thus the fraction of the true signal which 
we expect to measure with the OCRA-p system. We can 
then estimate the SZ signal which would be measured by 
an 'ideal' instrument with narrow beams at infinite separa- 
tion. (We note that the spherical model is non-ideal, but 
argue that in this context the adverse effects are small given 
that our observations are restricted to the central regions 
of the clusters. For more detailed analysis in the future we 
will adopt more realistic models, but do not feel that it is 
necessary here). We also account for the parallactic angle 
range for each set of observations. The results are recorded 
in Table It is interesting to note that F is very similar 
(and close to 0.5) for all clusters in the sample despite the 
large range of redshifts involved. This flatness of response is 
a generic feature of many practical SZ o bserving techniques 
(see, e.g.. iBirkinshaw fc Lancasterll2005l '). 



they offer the possibility of Fourier filtering. Techniques for 
implem enti ng this method are discussed in iGrainger et al.l 
(2002) and R cose et al.l (2002). For an instrument such as 
OCRA-p which operates a position switching strategy, radio 
sources can affect the data by producing a positive signal 
when they lie close to the cluster centre, or a rogue negative 
signal when they lie in the reference beam of the telescope. If 
we know the source flux densities at our observing frequency, 
we can correct the data and establish the true magnitude of 
the SZ effect for each cluster. 

No high-frequency survey of the radio sky exists, but 
for the four well-known clusters used for our commission- 
ing program we turn to the extensive SZ work reported 
in the literature. We are fortunate that the OVRO/BIMA 
group have observed all clusters in our sample at 28.5 GHz 
or 30 GHz, or both. We make use of the resu lts reported 



ICoorav et al l l| 19981 ) . iReese et all i|2002l ) and lCoble et~ail 
l|2006h . The primary beams of OVRO and BIMA are 4'. 2 
and 6'. 6 respectively, and so cover the full fields of rel- 
evance for OCRA-p . We also consult the all-sky NVSS 
ilCondon et all Il998l . 1.4 GHz) and GB6 l|Gregorv et all 
1 19961 . 4.85 GHz) catalogues, plu s the 5 GHz observations 
made by iMoffet fc Birkinshawl (|l989T) . We consider all 
sources within 5' of the pointing centre for each cluster. 

All potentially problematic sources are listed in Table 
[3] Additionally, the source positions relative to the tele- 
scope beams are plotted on RA-dec axes in Figure [1] The 
greyscale illustrates the arcs generated by the moving blank- 
sky beams, and the density of points is proportional to the 
observing time spent at each parallactic angle. The central 
circle marks the position of the 'on' beam at all times. 

We apply corrections to our SZ data based on the po- 
sitions and flux densities of the radio sources in each field. 
Where we have a measurement of a source at 30 GHz, we 
use the value directly. For the sources for which information 
is available at two or more lower frequencies, we fit a sim- 
ple spectrum and extrapolate to estimate the 30 GHz flux 
density. Although this method is non-ideal, we note that 
the error introduced by extrapolating from 28.5 to 30 GHz 
is likely to be very small. Where only one measurement is 
available, we adopt a cautious approach and run our cor- 
rection software twice: once using only the sources detected 
by OVRO/BIMA, and again using an upper limit on the 
30 GHz flux. Contamination through sources at 30 GHz is 
generally low for bright clusters, producing only small cor- 
rections to the SZ data, and thus introducing little bias as a 
result of the limitations of the methods employed. We do not 
make any allowance for source variability - this is discussed 
further in Section HI] 

We now examine the radio environment for each cluster 
in turn, and outline the procedure adopted for each. 



5.1 CI 0016+16 

We find no sources within 5' of CI 0016+16 and thus no 
correction to the SZ measurement is required. 



5 RADIO SOURCE CONTAMINATION 

Radio source contamination remains a significant problem 
for Sunyaev Zel'dovich observations. Here interferometers 
have a distinct advantage over single-dish experiments as 



5.2 MS 0451.6-0305 

For MS 0451.6-0305, source 1 lies inside the reference arc, 
and source 2 approximately one beamwidth outside. We are 
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Cluster 


r c 
(arcsec) 


P 


Offset 
(arcsec) 


F 


Refs 


CI 0016+16 


36.6±1.1 


0.697±0.010 


30.9 


0.52 


a 


MS 0451.6-0305 


45.1±1.2 


0.88±0.02 


17.6 


0.53 


b 


MS 1054.4-0321 


4+22.2 




21.7 


0.55 


c 


A2218 


56.0±0.84 


0.591±0.004 


28.4 


0.43 


b 



Table 2. /3-model parameters, estimated offset of pointing centre from apparent X-ray centre, and the e stima ted f raction of SZ flux 
measu red by the OCRA-p two-beam system. /3-model parameters taken from (a) Worrall & Birkinshaw] (l2003f) . fb) |Pe Filippis et al.l 
1120051 ). (c) iNeumann fc Arnaudl ll2000h 



fortunate to have a 30 GHz flux density measurement for 
source 1 which we use directly. Source 2 is detected only in 
the NVSS catalogue. It is conceivable that this source is lost 
in the noise of the OVRO/BIMA maps, and thus we run 
our correction software twice: once without the additional 
source, and again using an upper limit for its flux density. 
The value for the upper limit is taken as three times the RMS 
quoted for source 1, 0.78 ± 0.26 mjy. Our measurements are 
unaffected by the presence (or absence) of this source, since 
it lies outside the reference arc and has a low flux density. 

5.3 MS 1054.4-0321 

We find four possible contaminant sources in the MS 1054.4- 
0321 field. Source 1 lies close to the pointing centre. Source 

3 is found inside the reference arc so may also be a prob- 
lem, although it will enter with low weight. Sources 2 and 

4 lie outside and inside the reference arcs respectively. We 
have 28.5 GHz measurements for sources 1, 3, and 4 which 
we combine with the NVSS values at 1.4 GHz in order to 
estimate 30 GHz fluxes as before. For the remaining source, 
we apply a similar approach to that used for MS 0451.6- 
0305. We place an upper limit of 0.12 ± 0.04 mjy and re-run 
the source correction procedure. Again we find that the pres- 
ence/absence of this source does not affect our measurement, 
despite its position in the reference arc, presumably to the 
low value of the upper limit placed on its flux density. 

5.4 A2218 

A2218 has been observed at 30 GHz, and three sources are 
detected. We use these flux densities directly. 



6 RESULTS AND DISCUSSION 

Our SZ measurements are presented in Table 2] Column 
two contains the double-differenced values measured for the 
TARGET field for each cluster. Column three lists the same 
values for the TRAIL fields, which we note are generally con- 
sistent with zero. As mentioned in Section \4. II we feel that 
this eliminates the need for our third level of differencing and 
subsequently will use the TRAIL measurements as a check 
only. We illustrate our data via the cumulative histograms 
presented in Figure [2] For each cluster, the TARGET data 
are clearly shifted towards negative values relative to the 
TRAIL data, indicating clear detections. We apply a stan- 
dard two-sided Kolmogorov-Smirnov test and verify that 
at the 1% level, the measurements from the TARGET fields 
have different distributions from the TRAIL fields. 

The double-differenced values in column two of Table [4] 



represent only a fraction of the actual SZ signal (as discussed 
in Section 14. ip . We list flux density values after correction 
(see Section 14.41 for the instrumental signature and radio 
source contamination (see Section in column four. Cen- 
tral temperature decrements are summarised in column five, 
with the amount of telescope time listed in column six (note 
that this is roughly twice the actual integration time). We 
achieve clear detections of our four well-known clusters. 

To facilitate comparison, we also list our results along- 
side those quoted in Section [3] in Table [5] Our results 
for CI 0016+16, MS 0451.6-0305, and MS 1054.4-0321 are in 
good agreement with previous work, affirming that our ob- 
serving strategy, radio source treatment, and modeling are 
all effective. The one exception is A2218 which is rather high, 
although only ~ 1.5cr away from the other measurements. 
However, as discussed in Section ^. 41 this cluster displays ir- 
regular morphology. With the OCRA beam we may be mea- 
suring the SZ signal from a particularly high pressure region. 
Due to the extended nature of the SZ effect and ambiguities 
in defining cluster centres, it is likely that our pointing cen- 
tre differs from those used by the comparison experiments. 
Additionally, we note that we have not incorporated the pos- 
sibility of source variability in our analysis. We deduce that 
this is not a major factor here, given the good agreement 
between our measurements and those reported in the liter- 
ature, although this could perhaps provide an alternative 
explanation for the discrepancy noted above. 

The short integration time required to obtain these clear 
detections demonstrates the power of OCRA-p in SZ stud- 
ies. We stress that the values quoted in table 4 correspond 
to the total time required to take the data. The OCRA 
control system is rather inefficient (although recently im- 
proved) and thus the telescope spent more than half of this 
time off source. The predicted noise on a one second inte- 
gration of ~ 3.6 mjy degrades by a factor of ~ 2 to the mea- 
sured value of 7.0 mjy. Integrating over the full 'on-source' 
time of around 3 hours (8478 samples), we would expect 
this value to reduce to ~ 0.08 mjy, which is better than the 
measured ~ 0.4 mjy by a factor of 5. This degradation is 
likely due to non-optimal removal of atmospheric effects, 
which may partly be attributed to the poor telescope con- 
trol. With improvements to both the OCRA receiver and 
the telescope control system, we expect OCRA's sensitiv- 
ity to increase substantially in forthcoming observing runs. 
Over the current season, we have scheduled observations of 
a large, statistically complete, X-ray selected sample of clus- 
ters at z > 0.2. These observations will form an SZ catalogue 
with uniquely well-understood selection effect s. This will en - 
able us to test SZ/X-ray scaling relations (e.g. Cooray 1999). 
The SZ measurements will be combined with high-quality 
X-ray data in order to further test cluster models. We also 
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Cluster 


Source 


RA 


DEC 


S30.0 


•S28.5 


S3 


•Si. 4 


References 




No. 


(J2000) 


(J2000) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 




C10016+16 - - 


MS 0451.6-0305 


1 


04 54 22.1 


-03 01 25 


1.41 ±0.26 


1.86 ±0.26 




14.9 ±0.7 






2 


04 54 30.53 


-03 00 44.5 








14.6 ±0.6 


-, -, -, c 


MS 1054.4-0321 


1 


10 56 59.6 


-03 37 30 




0.90 ± 0.04 




14.1 ±0.9 


-j d, -, c 




2 


10 56 48.24 


-03 40 04.6 








4.6 ±0.4 


-, -, -, c 




3 


10 56 48.74 


-03 37 25.4 




1.4 ±0.04 




18.2 ± 1.0 


-, d, -, c 




4 


10 56 57.58 


-03 38 54.7 




0.4 ±0.04 




3.1 ±0.4 


-, d, -, c 


A2218 


1 


16 35 22.1 


66 13 23 


4.29 ±0.21 


4.43 ± 0.20 


2.8 ±0.2 


1.25 ±0.21 


a, a, b, b 




2 


16 35 47.7 


66 14 46 


1.36 ±0.10 


1.59 ± 0.11 


3.7 ±0.3 


18.0 ± 1.8 


a, a, b, c 




3 


16 36 16.0 


66 14 23 


2.41 ±0.29 


3.13 ± 0.30 


4.2 ±0.1 


13.3 ±0.6 


a, a, b, c 



Table 3. Radio sources as socia ted with each of the six clusters. Fluxes ta ken from (a) iReese e t alj d200 2fl. (b ) 
iMoffet fc Birkinshawl Jl989T) . (c) ICondon et ail dl998h . (d) ICoble et al.l J2006T) , (e) 4.85 GHz. lGregorv et alj (119991 . 
The source index refers to Figure [T] 



0016+16 



0451 




- c - 



x — offset/ arcmin 
(a) C10016±16 



x — offset/arcmin 
(b) MS 0451.6-0305 



1054 



A221S 








2 



x — offset/arcmin 
(c) MS 1054.4-0321 



-2 2 

x— offset/arcmin 
(d) A2218 



Figure 1. Radio sources in each cluster, shown relative to beams A and B, where the 'off' beam B traces arcs as the telescope tracks 
the cluster. The source numbers correspond to the details presented in Table [3] 
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Cluster TARGET flux TRAIL flux S sz AT Telescope Time 

(mjy) (mjy) QJy) (fiK) (Hours) 

CI 0016+16 -1.884 ±0.345 0.540 ± 0.411 -4.28 ±0.78 -1647 ± 302 11.5 
MS 0451.6-0305 -2.500 ± 0.425 0.446 ± 0.347 -4.05 ± 0.80 -1558 ± 309 10.5 
MS 1054.4-0321 -1.953 ± 0.406 0.331 ± 0.384 -4.48 ±0.74 -1722 ± 283 13.5 
A2218 -1.696 ±0.317 0.036 ± 0.286 -3.01 ±0.75 -1159 ± 288 1L0 

Table 4. Measured fluxes in the TARGET and TRAIL fields, difference and corrected measurements, and central temperature decrements. 




Figure 2. Cumulative histograms illustrating the raw data obtained through observing the Target and Trail fields for each cluster. Note 
the negative shift of the Target curve relative to the Trail curve, indicating clear SZ detections. In each case, the distributions are clearly 
separated by a standard Kolmogorov-Smirnov test. 



anticipate performing observations of high-redshift cluster 
candidates found by the XMM-Newt on Large Scale Struc- 
ture survey such as that reported bv iBremer et al.l |2006). 
The two samples will be at significantly different mean red- 
shifts, thus facilitating studies of the evolution of cluster 
properties. The OCRA-p commissioning observations also 
pave the way for SZ imaging. It is possible, by means of 
numerous pointings, to map clusters using OCRA-p. How- 
ever, a multiple-beam instrument such as OCRA-F would 
be far more efficient. The installation of OCRA-F on the 
telescope will allow us to make the first well-resolved maps 
of SZ clusters with OCRA-F after summer 2007. 

OCRA-F, and indeed the full, 100-beam instrument, 
will also be used to perform blind surveys for galaxy clusters. 
Some source correction will be possible via follow-up obser- 
vations of suspected sources present in the reference beam 
traces. Deeper observations, coupled with reference arc anal- 
ysis, would, allow contaminant sources to be detected and 



removed provided that the data are taken over a sufficiently 
large range of parallactic angles. Naturally, there may be 
some cases where sufficient additional information about 
possible contaminant sources is not available. For clusters as 
massive as the four discussed here, we do not expect radio 
sources to preclude detection. If we were to ignore the source 
corrections, all four clusters would still be unambiguously 
detected, and the temperature decrements would remain 
close to the values reported in the literature. C10016±16 has 
no source contamination so is unchanged. The AT values for 
MS 0451.6-0305 and A2218 become rather high as a result of 
switching onto sources in the reference arcs ( — 1814 ±308 /iK 
and -1517 ± 284 respectively) while MS 1054.4-0321 is 
low due to the central radio source ( — 1366 ± 284 /iK). 

As we move towards observing fainter clusters, we 
recognise that contamination by primordial CMB features 
could potentially be a problem, as noted in other work (e.g. 
lLancaster et al.ll2005l ). However, OCRA is a relatively high- 
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resolution instrument, probing arcminute scales which cor- 
respond to 3500 < i < 10000. At such high multipoles, for a 
concordance cosmology and including the effects of lensing, 
we expect the amplitude of the temperature power spectrum 
to be of the order of a few /iK. This would be well below the 
OCRA noise in surveys of scale ~ 100 square degrees with 
planned sensitivity of ~ 100 ^iK. 
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